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1. Introduction 

In the type of statistical models considered here, a fixed underlying graph carries 
a random realization of a structure, consisting of values (boolean, in a finite 
alphabet or scalar) carried by bonds (edges) or sites (vertices), and usually 
interacting via local rules. We will focus on the Ising and Potts models (and 
more generally cyclic and abelian spin models), where each site carries a random 
"spin" ; the discrete Gaussian free field, where a scalar height function is defined 
on sites; the random cluster model, where realizations are random subgraphs; 
and the six-vertex model, where each edge has a random orientation. 

The key question is to understand the large scale random structures generated 
by these local interactions. For instance, one may study the correlations between 
local configurations observed at large distance on the underlying graph. Wc will 
generically refer to (functions of) these local inputs as order variables. 

Historically, a fundamental tool, introduced by Kramers and Wannicr in 
the context of the Ising model, is duality, which maps the model defined on 
a (weighted) graph to the same model on another (weighted) graph. Of partic- 
ular interest are the "fixed points" of these duality transformations; they often 
coincide with a critical point of the model, at which one can observe large scale 
random structures. This type of duality is observed in a variety of models, and 
for abelian models can be phrased in terms of the Fourier- Pontryagin transform 
and the Poisson summation formula. 

Again in the context of the Ising model, KadanofF and Ccva showed that under 
duality, order variables are mapped to "disorder" variables, which represent a 
modification of the state space by the introduction of local "defects" ; moreover, 
the field-theoretic concept of fcrmion finds a combinatorial incarnation as the 



•Partially supported by NSF grant DMS-1005749 and the Alfred P. Sloan Foundation. 





imsart-ps ver. 2011/11/15 file: Cornell. tex date: December 20, 2011 



J. Dub edat/ Topics on abelian spin models and related problems 



1 



combination of an order and a disorder variable at microscopic distance. Such 
combinations are referred to as parafermionic or spinor variables. 

In these notes, we will describe duality, order, disorder, and spinor variables 
for abelian spin models and the discrete free field. We then discuss combinatorial 
mappings between the Potts model, the random cluster model, the six-vertex 
model, and the dimer model. 

For focus and (relative) sclf-containcdncss, we do not discuss here a number 
of major contiguous subjects. These include: phase transition, transfer matri- 
ces and Bethe ansatz ([2, 38]), Yang-Baxter solvability ([2]), infinite volume 
measures ([15, 18]) and limit shapes ([3, 27, 38]), scaling limits and Schramm- 
Locwner Evolution ([47, 43]). 

2. Spin models 
2.1. Ising model 
Definition 

In the Ising model ([32, 36, 19]) , a configuration on an underlying finite 
graph r = {y,E) consists of an assignment of a spin value a{v) G {1,-1} to 
each vertex v . The weight of a configuration is: 

w{(t) = exp(/3JeCT('(;)cr(w')) 

e—{vv')^E 

where Je 7^ is a coupling constant attached to the edge e G E and /3 > is 
the inverse temperature. Plainly the model (up to normalization) depends only 
on the edge weights: 

w{e) — exp(— 2/3 Je) 

It is occasionally convenient to allow w{e) = 0, i.e. Je = +00. Note that these 
weights are invariant under spin flip: w{—a) = w{a). One may also consider 
± boundary conditions, that is, fixing the values of spins to ±1 on some pre- 
scribed subsets of V (the configuration space then consists of spin collections 
{a{v))v£V s.t. (j{v) ~ ±1 for v G V^). One may also consider some segments 
on the boundary to be wired, viz. connected by edges with zero weight, so that 
the only contributing configurations are constant on these segments (a wired 
arc may be equivalently represented as an extended vertex). 
The partition function is 

Z = Z(r, ( Je)e6£;) = ^ wia) 

cr:V^{±l} 

We denote the expectation under the probability measure given by ¥{a ~ ctq) = 
w{ao)/Z by (•). 

A classical object of interest is the spin correlation 

{<t{vi) ■ ■ ■ a{vn)) 
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for vi, . . . ,Vn graph vertices (in particular when they are at large distance) and 
the spin variables a{v)''s are fundamental examples of order variables. 

The spin correlations depend on the boundary conditions. A wired arc may 
be seen as a boundary arc along which couplings are infinite (vanishing edge 
weights); in turn, ± boundary arcs may be represented in terms of wired arcs 
and spin variables. For example, if 7^,7" are two disjoint boundary arcs, we 
have 



where X is a generic random variable, the expectation on the LHS is for ± 
boundary conditions on 7*, and the expectations on the RHS are for 7+. 7" 
(separately) wired arcs. This shows that, allowing for general couplings, it is 
enough to treat free boundary conditions. 

From now on we will consider the planar case; Ft = {V^E"^) denotes the 
planar dual of F, so that faces of F correspond to vertices of F^^, and vice- versa 
(Figure 1). 



More precisely, it is rather convenient to think of vertices on a wired arc 
of the boundary as a single extended vertex; and to connect vertices on a free 
boundary or extended vertices to a vertex at infinity by edges with weight 1 
(i.e. with no interaction), li e £ we denote the corresponding edge of F^; 
for oriented edges, crosses e from right to left. 

Graphical expansions 

In the low-temperature regime (/? ^ 1), disagreements between neighboring 
spins are severely penalized. The low temperature expansion of the Ising model 
consists in mapping a spin configuration (cr„)„gv to a subgraph P = {V\Ep) 
of F^', where Ep = {e = [vv')^ E E^ : <T{v)a{v') = —1}- At low temperature 
this graph is typically sparse. The spins (ay) may be recovered from P up to 





Fig 1. Left: a portion of a planar graph F. Right: its dual F^ 



imsart-ps ver. 2011/11/15 file: cornell.tex date: December 20, 2011 



J. Dub edat/ Topics on abelian spin models and related problems 



3 



an overall spin flip (in the ease of free boundary eonditions). It is easily eheeked 
that P is an even degree subgraph of P; such subgraphs are called polygons. 

The (unnormalized) Ising measure on spins [a^) projects to a measure on 
polygons: 

w{P) = 1{P admissible} W w{e^) 
eeEp 

omitting a multiplicative factor 2 (for overall spin flip). 

Let us turn to the high temperature expansion. Since a{v)(j{v') = ±1, we may 
write 

exp(/3JeCr(v)cr(u')) = eosh(/3Je) (1 + tanh(/3 Je)cr(w)(T(u')) 
Let us define a dual edge weight: 

u>'(e) = tanh(/3Je) 

In the high temperature regime (/3 ^ 1), these dual edge weights are close to 1. 
Starting from the partition function, we obtain 

Z(r,(Je)e)= ^ W exp(-2^Jel{^(^,)^(„,)^_i}) 

<j:V^{±1} e={vv')eE 

= ncosh(/3Je) ^ n a + ^'ie)a{v)a{v')) 

eeE cr:V^{±l} e={vv')eE 

= ncosh(/3j,) En^'(^W"M^') 

eeE a:V~¥{±l} E„CE eeEo 

by fulling expanding the product. Then exchanging summations leads to 

z(r, iw,)e) = 2i^i n cosh(/3j,) n ^'(^) 

eG-E P polygon eei?p 

Then comparing high and low temperature expansions, we have the Kramers- 
Wannier duality ([29, 30]) for partition functions: 

Z(r, (Je)e) = 2l^l ( n COSh(/3Je) J Z{T\ (4)) 

where 

= tanh(/3Je) 

or w' {e\ — T-r^^T^. Note that w i-^ is involutive. It exchanges and 1, 
SO that wired boundary components are exchanged with free boundary compo- 
nents. To retain positive weights after the duality transformation, we need to 
start from ferromagnetic couplings, i.e. Je > 0. 

If F is (a portion of) the square lattice, is identical to F, up to boundary 
modifications. For instance, one can consider F to be a rectangle with half-wired 
and half-free boundary eonditions, in which case F^ is exactly isomorphic (as a 
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weighted graph) to F. Thus for the square lattice the fixed point Wgd = V2 — 1 
of the duality mapping w i— > is self-dual, in the sense that it is fixed for 
Kramers- Wannier duality (at least "in the bulk" ) . The low-temperature regime 
w < Wsd and the high-temperature regime w > Wsd are exchanged by duality. 

Disorder variables. 

Let us now consider the effect of Kramers- Wannier on spin correlations 
{a{vi) ■ ■ ■ a{vn)) ■ In the high temperature expansion (and periodic boundary 
conditions for simplicity), we can write: 

Z(r,(J,),)(c7(«i)...a(z;„)) = n e'^^=^('')'^(''') 

and then 

Z(r, (Je)e) ,^ , \^ ^ 

n,^cosh(/3Je)^^"^ • • • -.J = 2^ 2^ 

= 21^1 E n -'(^) 

Per eeEp 

where the sum bears on subgraphs P C F which have odd degree at {vi, . . . , «„} 
and even degree elsewhere (by a <-> —a invariance, the correlator is zero if n 
is odd, in which case the sum is empty). For periodic boundary conditions 
(i.e. when F is embedded on a torus), one requires additionally that any non- 
contractible (with respect to the torus) path on F^ crosses an even number of 
edges of P. 

Such a polygon with defects at {vi, ...,«„} does not correspond to a spin 
configuration on F^^. Following Kadanoff and Ceva ([23]), this motivates the 
introduction of disorder variables, which encode a modification of the state 
space. 

The data of {av)vev is equivalent, up to global spin flip, to the data {da{e))eeE, 
where da{vv') = (t{v')(t{v)~^ . li w : E ^ we deflne: 

dc^if) = H uie) 
eedf 

where / is a face of F and df is its (counterclockwise oriented) boundary. Up to 
spin flip, spin conflgurations (tT^)i,gy correspond to closed currents: {(we)eG-E : 
doj = 1}. (In the periodic case, there is an additional condition: nee7'^(^) ~ ^ 
if 7 is a non-contractible path). Plainly the Gibbs weights can be written as 
functions of these currents. 

Introducing disorder variables /^(/i), ....... /i(/n) (/i, . . . , /„ faces of F) con- 
sists in modifying the state space to: 
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with weights 

w{uj) = w{e) 

We assume that n is even (otherwise the state space is empty, as O/grt '^if) ~ 
1 for all w). If u, v' G V, one can define a{v')a{v) = Yl^^^ , a;(e) where jv^v' 
is a path from v to w' on F. Notice however that this definition depends on 
the choice of 7. Indeed, if 7' is another such path, 0667 '^(^) ~ ^nee7'-^(^) 
depending on the parity of the number of fi's enclosed by the loop 7'7~^. 

One can trivialize this data by choosing ^ "defect lines" (paths on T^) joining 
the /i's pairwisc. Then one can find a : V {±1} such that uj{e) = (icr(e) if 
e does not cross one the defect lines. In terms of this trivialization, wc have an 
Ising model with antiferromagnetic (negative) couplings — Je for e G -B crossing 
a defect line. 

We can now define mixed order-disorder correlators: 

i j u>;E->{±l} e67 e£E:ui(e)=-l 

where {vi, . . . , V2m} G T, F = {/i, . . . , /2n} C F^^, 7 is a union of m paths on F 
with endpoints vi, . . . , V2m- The sign of the correlator depends on this implicit 
choice of paths (or alternatively of defect lines, in which case one simply requires 
that 7 does not intersect defect lines). 

In this case the Kramers- Wannier duality ([29, 30, 23]) yields: 

i 3 i 3 

where the couplings on Ft are in duality with those on F. 

Kadanoff and Ceva ([23]) identified the combination of an order variable a 
and a disorder variable /i at microscopic distance as a discrete version of the 
field-theoretic notion oi fermion. In particular one denotes tpyf = a^jif where 
i; G y is on the boundary of the face /. 

2.2. Abelian spin models 

The Ising model may be generalized in several directions; the first we consider 
here is that of abelian spin models, in which the spin variable at each vertex 
takes values in a fixed finite abelian group G. In the case G = {±1}, one recovers 
the Ising model. 

Given a graph F = [V, E) and a finite abelian group G, a configuration is a 
mapping a :V ^ G. For each e £ E,we have a weight function We ■ G ^ [0, 00), 
which is symmetric, i.e. We{g~^) = We{g) (relaxing this condition leads to chiral 
models and requires the underlying graph F to be oriented). Then 

w{a) = 'w{(j{v')a{v)~^) 

e—{vv')^E 
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is the weight of a general configuration. For free boundary conditions, we have 
'w{pu) = w(a) for any p <E G. 

If G = TLjqL, we have a clock model. If furthermore We = a + 61{i)., we get the 
q-state Potts model. If g = 4 and We is a general symmetric weight, we obtain 
the Ashkin-Teller model ([2], as phrased by Fan). 

Let G = Hom(G, C*) be the dual of G [G and G are isomorphic, though 
not canonically) . The basic order variables are of type x(ct(z;)), x € G, with 
associated correlators: 

{\{x^{cr{v^))) 

i 

where Vi £ V, Xi G G. For free boundary conditions, this is zero unless Yii Xi = 
1. 

Kramers- Wannier duality has been extended successively to various models, 
see in particular [48, 41]. By Fourier- Pontryagin duality ([40]), one can write: 

where We is the Fourier-Pontryagin transform 

|G|i/2 

At the level of partition functions, we have the summation: 

z(r,G,K)e)= E n M^W)a{v)-') 

a:V^G e={vv')eE 
= |G|-^I^I E n E ^e{Xe)XeHv')a{v)-') 

= |G|-^I^I E n We{Xe)XeHv')aivr') 
a-.v-^G e—(vv') 

= \Gr\-m E n E^(xe) 

X■■E-^G■.d*x=l e={vv') cr„(EG 

where d * x(^) — Ylv'~v Xe- Let us identify x '■ E ^ G with x ■ ^ G via 
x(e^) = x(e) for e £ E. Then d * x = 1 is equivalent (in the simply connected 
case) to the existence ofa -.V^ ^ G such that xiff) = 9{f')9{f)~^ ■ (Implicitly, 
we fixed a reference orientation for edges of F^; by symmetry of the weight, one 
can write each term We{(T{v')a{v)~^) so that v (resp. v') is the left (resp. right) 
vertex wrt the oriented edge of F^^ separating them). We conclude that: 

Z(F, G, (w,),) = |G|l^l-^l^l+iZ(Ft, G, {we)e) 

Disorder variables are indexed by a face f € and an element g G G, and 
denoted Pg{f). The corresponding modified state space consists of G- valued 1- 
forms w : £ — > G defined on oriented edges and antisymmetric in the sense 
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that ojixij) — uj{yi) ^. Wc define doj{f) — H-g-Ga/ "^C^)' where df is the coun- 
terclockwise oriented boundary of /. Then the state space corresponding to a 
disorder niMsi(/») is: 

{w : du = Y[9^^'} 

i 

with weight w{u}) = '[l^_,^^We{oj{e)) (there is an overall factor \G\ compared 
with the earlier normalization) . Again this may be trivialized by fixing a defect 
line going through all /i's and modifying couplings along this line. Remark that 
Y[i fi^^ (otherwise the state space is empty) . For simplicity let us consider a 
pair , i<^ivi) of spin variables, and a pair /Xg(/2), /Xg~i (/i) of disorder 

variables. Let us fix non intersecting paths 7 (resp. 7''') on F (resp. F^) from vi 
to V2 (resp. from /i to /2). We can write a mixed correlator: 

Z{x{a{v2))x-\<7{vi))fig{h)fig-.{fi)) - 

E n xi'yiv'Mvy') n w',iaiv')aiv)-') 

where w'^ = if e does not cross 7'^ and We(-) = We{g-) otherwise (in which 
case we orient e from right to left of 7^'). Repeating the argument above leads 
to 

(Xo(cr(w2))X(7^(cr(wi))/ig(/2)Aig-i(/i))r = (Mxo("2)Mx-i(^'i)x/2(3)x/i(3"^))rt 

where the x/'s are the dual order variables; weights on Ft are Fourier coefficients 
of those on F. This extends to any number of insertions. 

By analogy with the Ising model, on the square lattice we are particularly 
interested in self-dual weights, viz. weights such that {? cx w o </> for some iso- 
morphism : G — ^ G. If we assume weights are nonnegative, by Parseval we 
have w ~ w o (j). 

Cyclic models. 

Let us specialize to the cyclic case: G — Z/qZ. Vq — {z G C : z'^ = 1}. Let 
.^0 = exp(2i7r/(7). Identifying G with G in the usual way, the Fourier transform 
is written: 

Finding a self-dual weight consists in finding eigenvectors for J-". This operator 
is the discrete Fourier transform (DFT), involved in fast Fourier transform al- 
gorithm. The question of its diagonalization is classical, if rather delicate. By 
Fourier inversion, the Fourier transform J" : — C*^ is of order 4: .F^ = Id, 
and consequently the eigenvalues of J-' are {±l,±i}. It is known that the mul- 
tiplicity of the eigenvalue 1 is [|J -1-1. 

First we look for self-dual Potts model weights. In this case: w cx 1 -I- qSq. 
Plainly = -j^ + ^5q. The fixed point equation Tw = w is thus solved 
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for a — yjq. For Ising (q ~ 2), we recover the self-dual weight function (up to 
multiplicative constant). 

As soon as g > 4, the multiplicity of 1 is > 2. Similarly to [20], one can look 
for a nice operator D G End(C'^) (almost) commuting with F. An educated 
guess gives: 

D = {a + bx'^)R+ic + dx'^) 

where R is the right shift: (i?/)(fc) = f{k + 1), and x € C-^ is the character: 
x(fc) = ^0- Classically, J^R = x^ (x acting by pointwise multiplication) and 
TX^^ = R^- Besides Rx = ^oxR- Then 

J'D = {ax + b^oxR + C + dR)T 

Consequently, J^D = XxDJ" (A G C) iff: 

(a, b, c, d) ~ Xib^Oi d, a, c) 

that is: = ^q^, {a,b,c,d) oc (1, A'^, A, A^). An element w of KerZ? satisfies a 
trivial recursion: 

wik + 1)(1 + A^^o"^-) + wik){X + A'Co"^') = 
To ensure positivity, we choose A = — exp(— z7r/2(7); then we find 

sin(2^ + f ) 

and check that w(g) — u'(O) ~ 1. Then KerZ? is a line which is invariant under 
J^, and thus w = is an invariant weight. This is the Fateev-Zamolodchikov 
point ([11]) for the Z/gZ model (identified in [11] as part of a family of weights 
satisfying the Yang-Baxter relations). 

Ashkin- Teller model. 

In the Ashkin- Teller model ([2], as phrased by Fan), each vertex v G V carries 
a pair of spins (cr^, py) G {±1}^; the interaction weight on the edge e = (vv') is: 

exp(/3(JeO-t,(T„/ + J'^PvPy' + J'ja^ay'pypy')) 

We then recognize the abelian spin model with G ~ {±1}^ (the simplest non 
cyclic finite abelian group). Identifying G with G we get (here (£1,62) G G) 

(J-/)(ei, £2) = I (/(1, 1) + £i/(-l, 1) + £2/(1, -1) + ei£2/(-l, -1)) 

leading to the following simple self-duality condition for the weight w : G — > M 
(i.e. satisfying w — J-w): 

w(l, 1) = wi-1, 1) + u>(-l, 1) + w(-l, -1) 
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The basic order variables are cr, p, ap and there are naturaUy corresponding dis- 
order variables p.^Vjpv. Remark that in the isotropic ease J = J' , the model 
becomes a Z/4Z model. 

Parafermions. 

Let us go back to the general case of an abelian model with spin group G. 
Order (resp. disorder) variables are indexed by G (resp. G). By analogy with 
the Ising model, one can consider a parafermion (or spinor) 

V'?f^° = xoKKo(/) 

where (xo, go) ^ G y- G (and are fixed and omitted from now on), w G F, / e F''' 
adjacent. When tracking an order variable Xo{iy) along a cycle around a disorder 
variable Mgo(/) (or vice versa), one picks a phase Xo(5o)^^- In particular, when 
the spinor rotates unto itself (eg fix / and follow v along df), it is multiplied by 
Xo(5o)^^- This may be formalized as follows. Consider a graph F' with vertices at 
midpoints of segments [w/], u G F, / G F^^, w on the boundary of / (denote such 
a vertex (w/)); edges of F' are such that {vf) ~ («'/) if v,v' are consecutive 
vertices on df and dually {vf) ^ {vf) if w G 9/ n df . The "line bundle" 
®(vf)^T' '^("/) is equipped with a connection (in the terminology of [25] eg), 
i.e. a collection of isomorphisms : Cu — > Cu' for adjacent vertices u,u' in 

F'. These are defined by: 

'i^(vf),(v'l){'A = cxp ^zsarg j _ ^{vf),{vf'){z) = cxp (^s&vg^^—j^ z 

where s is the spin, and arguments are chosen in (— tt, tt) (so that 4>u.u' ~ 4'u^u)- 
Let us consider a general correlator: 

n m 

F{vojo) = ii'vofo n ^'(^*) n (-^j)) 
i=i j=i 

We regard all insertions as fixed except u, /, which are adjacent. In order for 
this correlator to be non-trivial, we need Xorii>iXi — ^q-, 9oY[j>i9j — ^G- 
To assign an actual value to the correlator, we also have to fix a defect line (or 
defect tree) connecting all disorder operators on F^ . Alternatively, one can think 
of F as being multivalued (with a phase ambiguity) or as a section of the line 
bundle described above (with the necessary modifications around Vi,fj). 

Consider the local situation around an edge e = {vv') G E and its dual edge 
(//') = (ww')^ G i?^; this gives four possible locations for the parafermion. We 
assume that the defect line 7 chosen for / is the concatenation of (//') and the 
defect line 7' for /' (which has /' as an endpoint). 

For 5 G G, let us consider the partial (modified) partition function: 

2^g= ^<yiv)=lMy')=a]lx^{'^{v^)) Yl ^'dv^''^) 

(7:V~¥G i>l e=(xy)eE.e^{vv') 
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where w' is the weight function modified along 7' in order to encode the disorder 
Mso(/') Y{j>i l^gAfj)- Th'^n up to a factor \G\ 

F{v, /) = ^ We{9go)Zg F{v, /') = ^ We{g)Zg 
geG geG 

F{v', /) = E Xo{9)we{9go)Zg F{v', /') ^ ^ xo(ff)w^e(3)2<, 

see gGG 

Let us assume that the weight w = We satisfies the following equation on C*^: 

i{a-bxo)Oo-{c-dxo))w:^0 (1) 
where is the shift: {9ow){g) = w{ggo)- Then 

- aF{v, /) + bF{v', f) + cF{v, /') - dF{v\f) - 0, (2) 
the local parafermionic equation ([43, 44, 39, 21]). Equation 1 implies 

k 

where z = Xo(3o)- In particular, if r is the order of go 



If r' is the order of z = xo(.9o) in U, we have: nj=o^(^ ~ -^"') = ^ 1- Conse- 
quently we need: 



Symmetry of the weight ^(.g ) = w(.g) yields further conditions, such as: 



n*^ c — dz-' c — dz ^ ^ ^ 
a ~ bz3 ' a - 6z--?'-i ^ 



coming from wlg^) = w(.go (assuming weights are nonzero). Hence 



c — dz*^ 6z — az'' 



Two homographies that agree on three points arc equal. Thus if r' > 3, we have 
(a, &, c, c?) oc {dz~^,c,bz~^,a), ie {a,b,c,d) cx {uX^^,u^^X,u^^X~^,uX) where 
u € C*, = z. In order for weights to be positive, we take u unitary and 
recover the FZ weights ([11]): 

'^"l Qinr^ -I- -^^ 

= n -^^feiHr 

j=o snu^:^ — j 

with an additional "anisotropy" parameter 6 compared to the self-dual case. (In 
the cyclic case, arguing as before we find that Twg oc wiL-e)- 
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3. Discrete Gaussian Free Field 

We have considered abelian spin models with values in a finite abelian group. 
More generally one may consider a locally compact abelian group (the natural 
set-up for Fourier- Pontryagin duality). Let us consider scalar fields, i.e. fields 
taking values in R. 

As before, we consider a finite connected planar graph F — {y,E). It is 
somewhat convenient to designate a non-empty subset 9 of F as the boundary. 
We consider discrete scalar fields, i.e. functions : V — >■ M. The values of on 9 
are fixed to some prescribed 4'\s. This defines an affine state space. Let us also 
fix positive conductances (ce)eeB- 

The Discrete Gaussian Free Field (DGFF) on (F, (ce)) with Dirichlet bound- 
ary condition on d is the Gaussian variable with measure proportional to 

n exp(-|(0K)-</>H)^) n ''^^"^ 



e=(vv')eE v<£V\d ^ 

For / : F — > R, its Dirichlet energy can be written: 

by a discrete Green's formula argument, where ZV is the weighted graph (posi- 
tive) Laplacian: 

(^/)(«)= Ec(™,)(/W"/(^')) 

Under our assumptions, 

Ar:{f:V^RJig = 0}^{f:V^ M,/|a = 0} 

is indeed invertible and the Dirichlet energy is a positive definite form on {/ e 
R^, /|a = 0}. Consequently the above Gaussian measure is finite. The mean field 
(po is the one minimizing the Dirichlet energy in the affine space {/ S R.^, f\d = 
<pd}, i-e. the solution of the Dirichlet problem: (ZV/)y\a = 0, f^g = (j)Q. 

The covariance kernel C{v,v') = Cov{<j>{v)<j>{v')) is the Green kernel G for 
/V, viz. is characterized by A"G(., v') — S^' mV\d, G{v, v') = for w e 9. We 
deduce the characteristic function: 

E(exp(i ajcpivj))) = exp{i ^ ajcpoivj)) exp(-^ ^ ajakG{vj ,Vk)) 

3 3 j,k 

The local variables exp(za(/)(?;)) are called electric operators, where a is the 
charge; these are the main order variables. 

In order to introduce disorder variables, it is convenient to rephrase the prob- 
lem in terms of the current J{vv') = cf>(v') — 4>{v), a graph 1-form (in f2^(r)). 
Let us also assume that there is no boundary. Then the Dirichlet energy is well- 
defined in terms of J: \ ^^Ce{J{e))'^ . As a state space, we consider the space 
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of f -forms such that dJ ~ m where m = rrij 1 S MX is a fixed magnetic 
charge distribution {rrij is the magnetic charge positioned at the face fj); here 
dJ{f) = "^g^Qf Jie)- The Gaussian variable J induced on this af&ne space by 
the Dirichlet energy can be written J = + Jg where J/j, the mean, is the 
current of minimal energy with given magnetic charge distribution and Jq is a 
centered field with covariance kernel: 

Cov{J{viV2),J{v3V4)) = G{V2,V4) + G{vi,V3) - G(W2, W3) - G(wi, W4) 

(In the absence of a boundary, G{.,v) is well defined only modulo an additive 
constant, which is enough to define the RHS). 
li J = d(j) is an exact f-form, we have 

Ce JeWe = ^ (/)(«) ^ Cy^'Loivv') 

which shows that the orthorgonal in ^ CgJe) of exact I-form are har- 

monic 1-forms, i.e. those satisfying J^v'r^v Cvv'^^iw') = for all v. It is con- 
venient to define a discrete Hodge star operator * : ri^(r) — >■ ri^(r''') defined 
by {*oj){{vv')^) = Cvv'Uj(vv') ([33]). We regard as a weighted graph with 
Cet = (ce)^^, so that *^ = — Id and * is isometric. Clearly, * maps harmonic 
forms on F to closed forms on F^', and vice- versa. Thus J S n^{E) can be de- 
composed uniquely as J = d(f> + *d'ip where dcj) g ^^^(F) and dip g ri^(F^) closed; 
this Hodge decomposition is orthogonal in L'^{E,'^CeSe)- By analogy with the 
continuous case, we write (J, J') l^(e) = / J/\*J' (see [33]). Moreover, ZV = d*d 
(as usual identifying vertices of F with faces of F'''). Consequently, the form J/j 
with minimal Dirichlet energy among those such that dJ = m is: 

V/6Vt J 
where is the Green kernel on Ft. 
Duality. 

Let us consider abelian duality for the centered DGFF. At this point we 
could write mixed electric-magnetic correlators in terms of G, G^ and see duality 
appear in these explicit formulae. For ease of comparison with the finite case, 
we phrase duality in terms of Fourier transform. 

In general one can fix an electric charge distribution (a-u)v<£V and a magnetic 
charge distribution {mf)f^p, such that J^v'^v = 0, J2f''^f = 0- For simphcity, 
let us consider two pairs of charges: electric charges ±a at w±, and magnetic 
charges ±m at f±. If v± circles counterclockwise around /+, the correlator picks 
up a factor exp(±iam); assuming am € 27rZ removes multivaluedness issues. In 
the general case, let us fix paths 7,7^ on T,T^ with endpoints v±,f±. We are 
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interested in (cxp(iQ! J^J)), where (.) is the unnorniahzed expectation for the 
free field current J on F with mean Jm- We have the representation: 

i) A *dJl 

J-f J 

where = Et,ey w). 

It is technicaUy convenient to temporarily mollify the hard constraint dJ = m 
and dampen the electric correlator by considering the Gaussian integral: 

r f £ T-r e~^^"'°^^ 

I exp(iQ; \ d4>- -||d(/)|p)exp( — 1| * dtp - Jm.\?) TT —==rdJe 

2 2 ^^^^ V27rcr^ 



where £ is a small positive parameter, J = dcj) + *dip (||.|| is the weighted 
norm, for which the Hodge decomposition is orthogonal). By Fourier inversion, 
this is written as 



/ 



1 {'J )^ 

exp{~'—\\d4>~dJi\\^)cxpii [ d^AdJ„,~'-\\di^\\^)r\^^=^dl 
2 J 2 ^27rce 



where J = d(f> + *dtp is the Hodge decomposition on Ft. 

Letting e \ 0, we conclude that the discrete free fields on F,rt (with recip- 
rocal conductances) are in duality, in which magnetic and electric charges are 
exchanged. 



Compactification. 



In the case where the underlying graph F is embedded on a torus S = C/L, 
if we consider the Gaussian measure on closed currents J : dJ — induced by 
the Dirichlet energy, we gain two "topological" marginals: {Jj^J,J^J), where 
A,B are two standard cycles generating the homology of E. Given a,/?, there 
is a unique closed, harmonic form dwa,p {doJa^ = 0, c? * uja.p = 0) such that 
{Jj^Ua^ls, JgUJa,i3) = (a,/3). The space of 1-forms decomposes as an orthogonal 
sum of exact, coexact, and closed and harmonic forms: this is (a discrete version 
of) the Hodge decomposition. Thus an instance of a free field current J can be 
decomposed as J = dcj) + uj, where is a scalar field (on V, defined modulo 
additive constant) and w, the instanton component, is a closed harmonic form; 
these two summands are independent. 

It is classical ([14, 4]) to consider a compactified version of the field, in which 
in the periods (/^ J, /g J) are constrained to take values in a prescribed lat- 
tice 27rrZ (r is the compactification radius). This has no effect on the scalar 
component of the field; the instanton becomes a discrete variable supported on 
{w27rmr,27mr}, with wcights proportioual to the Dirichlet energy of these forms 
(a "discrete Gaussian" variable). 

If A is a lattice in a Euclidean space {V, (., .)), its dual lattice A' (y g A' iff 
{x,y) e Z for all x G A), we have the following version of the Poisson summation 
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formula: 

5^/(x)=Vol(F/A)-i5]/(y) 

in the normalization /(y) = e^^"^^'^^/(a;)(ia;. In particular for /(a:) = e^'^*"^" 
we get 

Let us specialize this to V the two-dimensional space of closed harmonic forms 
(which represents H^{Yi)) with norm given by 2it times the Dirichlet energy; 
A = {cj : Jj^gOj E 27rrZ}. In order to identify A', we introduce the closed 

harmonic 1-forms oj^ ^ on such that Jj^ ~ ^^^^ ^^'^ recall the bilinear 

relation: 

{0Ja,b,*0jlj = / OJa.bAi^lii^ f ha.bi^l ^ = ad - be 
JT, ' JdF 

where _F is a fundamental domain bounded by the cycles a, b and dha^b = ^a,b in 
F. We conclude *A' = {uj^ : ^u}^ £ 2Trr~^Z}. Thus if we write the instanton 
partition function 

weA 

which represents the (multiplicative, by independence) contribution of the in- 
stanton component to the partition of the compactified free field at radius r, we 
have 

up to an elementary multiplicative factor. Thus duality on compactified free 
fields has the effect of inverting the compactification radius ([14]). This is a 
simple example of T -duality. 



Gaussian free field. 



The continuous Gaussian (or massless) free field on, say, a torus S is the 
Gaussian field with action 

5(0) = I 

(g > the coupling constant)^ i.e. the centered Gaussian field with covariance 
kernel ^(— A)^^ where we consider (—A) as an invertible operator on zero-mean 
functions on S ([42, 16]). In the abstract Wiener space approach, one thinks of 
the free field as a random element of a large enough Banach space, typically a 
negative index Sobolev space. Plainly it may also be considered on the plane, 
finite domains, Euclidean spaces of other dimensions and Riemannian manifold. 
The compactified free field with values in M/27rrZ is defined as above as the 
sum of a scalar free field and an instanton component, which is supported on a 
lattice of harmonic 1-forms. 
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Let US briefly discuss magnetic and electric charges in the plane ([5], [14]). 
In this case, the covariance kernel is C(z,u;) = g^^\og\z — w\ (with a slight 
abuse of terminology, as we have to quotient by constant functions). The electric 
operators are written formally Oe{z) = exp(ie0(0)) (e is the electric charge); this 
is somewhat problematic, as a realization of the field (j) is only a distribution. In 
the standard regularization, one disregards the diverging self-energy and write: 

(:aWO-e(u.) :) = \z-w\-'^ 

where the colons denote the regularization procedure. Alternatively, one may 
replace (^{z) with its average on D{z, 5) and define 

(: a(2)0-eM :) = limr^ (exp(ze ( <t> ~ ^ I <IA) 

-5^" \ y^o^ Jd(z,s) Jd{w,s) )I 

In the presence of a magnetic charge m at z (represented by ©^(z)), the field 
increases by 2rm: when tracked along a counterclockwise around z. For a pair 
of insertions Om{z)0^m{w)., the average field is TO(arg(. — z) — arg(. — w)), 
the harmonic conjugate of m(log |. — z| — log |. — w\). Its Dirichlet energy has a 
logarithmic blow-up at singularities; a standard finite part regularization gives 
the expression 

{■Oe{z)0.e{w) :) = \z-w\-^"s 
Combining these two elements, we obtain: 

(: 0,{z2)0^e{zi)Ora{w2)0^,^{wi) :) = 

\^Zi — Z2\~^\wi — yj^\-™-^ a ^i™.e(a,):^{z2-W2)-a.J:s(^2-wi)-a.ig(zi-Wi)+a.rg(zi-wi)) 

To obtain "spinor" variables Oenn we coalesce electric and magnetic charges: 
''Oem{z) = limz^wOe{z)Omiwy\ Specifically take Zi = Wi + 5ui, \ui\ = 1, 
(5 \ 0. Then 

(: Oem{w2)0-e,-m[wi) :) = |u;2„it;^ | " ^ -"'Se-^"^" »2-«;i ) giem(arg(«i )-harg(«2) 

The last part is a manifestation of the spinor nature of these variables. Fixing 
a reference direction for mi,M2, this may be rewritten as: 



(: Oem(w2)C-e,-Tn(wi) :) = (W2 -U^l) ^9 2" {w2 ~ Wi) " 



In this context, one recovers the duality between (properly renormalized) 
electric and magnetic operators Oem ^ C2Tn,e/2 (and g o | for couplings), and 
pairs of reciprocal radii in the compactified case ([14, 4]). 

4. Random-cluster model 

Let us consider the Ising model on a graph F = iV^E). The weight associated 
to the edge e and the spin configuration {av)v^v is exp(/3Je(T„(T„'). As we have 
seen, writing this weight as 

exp(/3Je(Tt,fTi,') = cosh(/3Je) (1 + (Tt,CT„' tanh(/3Je)) 
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leads to the high-temperature graphical representation of the Ising model. It 
may also be written as 

exp(/3Jea,a,0 = e^^'^^ (l + (5,„,,^, [e^^^^ - 1)) 

which leads to the random- cluster representation ([18]), as we now explain. Set 
w{e) = e^^'^' — 1. The partition function reads: 

Z= cMPJe(Tv(Jv')=\{e-^^' {l + 5^^^^^,w{e)) 

a:V^{±l} eGV cr:y->{±l} 

= n^-''= E E n ^^^^^.'-(^) 

eGV (j:V^{±1} Eo^E e={vv')eEa 

We notice that Y.a:V^{±i}Il{vv')<^Eo^'yi.v)M(v-) = 2^'-^"\ where C{Eo) is the 
number of connected components of the subgraph Fq = {V,Eq). Hence 

EoCE eeEo 

More generally, one may consider a g-state Potts model with edge weight e'^^''"^"'"''' 
which leads to the following random cluster representation of the partition func- 
tion: 

^ = E 1''^''°^ n ^(^) 

-EoC-E eeEo 

with w{e) = e'^^'^^ — 1. Plainly, the RHS may be interpreted as the partition 
function of a model where the configuration space consists of subgraphs Fq = 
(y, Eq) of F and the configuration weight of Fq is proportional to 

eG-Eo 

This is the random- cluster model or Fortuin-Kasteleyn (FK) percolation ([18]), 
and is defined for any q > (whereas the Potts model is restricted to q £ W). 

One may also track natural correlators, such as spin correlators, through the 
random-cluster representation. We sec that in the q-state Potts model, 

(CT(w2)a-(wi)"^)potts = ^fk{vi O V2) 

where vi V2 means that vi , V2 belong to the same connected component 
(cluster) in the random subgraph induced by FK-percolation, and Potts spins 
arc regarded as elements of = {z G U : = 1}. By the &q symmetry of 
the Potts model, we have (5o-(„)4)potts = ^ and {a{v2)cr{vi)~'^)potts =P + jE^ 
where p = {Sa{vi),a(v2)); then 

Cov((5o.(t,i)4,(5o-(„2),i) = - — \ = ^-2-(o-(w2)o-(wi)~^)potts = ^— 2— Pf/<'(wi ^ v-i) 
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Higher order spin correlations are easiest to understand with the Edwards- 
Sokal couphng ([7]) between the g-Potts model and the corresponding FK per- 
colation. Starting from an FK configuration, one chooses a spin uniformly at 
random for each cluster and assign it to each vertex in this cluster. Then the 
resulting spin configuration is q-Potts distributed. Consequently, a general spin 
correlator Xj{^{''^j)))'Potts may be expressed in terms of the connectivity of 
wi , . . . , u„ in the FK percolation subgraph (which induces a random partition 
of {ui, . . . ,u„}). 

Duality. 

In order to discuss duality, we need to specialize to planar graphs. In perco- 
lation and related models ([17]), it is standard to introduce a dual configuration 
eI C E'^ defined by e^' G E'jJ iff e ^ So, so that |£'o| + \eI \ = \E\. Euler's formula 
for Fo {V,E(i) reads: 

\V\-\EQ\ + \Fa\ = l + C{E^) 

where Fq designates the faces of Fo (on the Riemann sphere) and is easily seen 
to be in bijection with the connected components of fJ = {V\ e\). Hence 

qCiEo) -Q ^(g) ^ L\V\^\E\^l -Q ^CiEl) -Q 9 

eG-Eo V eG_E / et(E_E^ 

Hence by defining a dual edge weight w{e') — —j-^, we see that Eq ^ e1 is a 
measure preserving correspondence between FK percolation on (F, {w{e))f,(£E) 
and (F, {w{e'^))^f^Ei)- 

In the case q £ {2,3,...}, this is consistent with the abelian duality for 
the Potts model, seen as a Z/qZ model. Indeed, if the edge weight function 



its Fourier transform is: 

yielding back w{e^) = 
Loop representation. 



We OC 1 -I- WeSl 



In turn, the random cluster representation (in the planar case) may be 
mapped to a dense loop (or fully-packed loop) representation (in the standard 
case of the square lattice with self-dual weights, this is the dense 0(n) model at 
n = y/q), which we now describe. Attached to the pair (F,F''') of graphs on the 
Riemann sphere, one constructs the quadrangulation <0> (the "diamond graph" , 
[33]) with vertices 1^ = V UV^ and edges E^^ defined by {vf) S is w is a 
vertex of V on the boundary oi f E (identified with faces of F). Then D(}, 
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the derived (or medial) graph of is specified by Vd<> ~ and (6162) S -Eo 
if 61,62 are two consecutive edges on a face of (}. (This is a planar version of 
the derived graph). Note that exchanging the roles of r,r^^ leads to the same 
0, DO. 

To a pair of dual configurations Eq d E, e\ d E\ one associates E'q C i?_DO, 
the set of all edges in DO that do not intersect edges in Eq,e\. Specifically, 
if {V1V2) e E, (/1/2) e E\ then either {V1V2) & Eq and thus ((/iwi), (/i'U2)), 
{{f2Vi), if2V2)) e E'^- or (/1/2) e El and thus ((/i«i), (/2«i)), {{hv2), (/2«2)) S 
i?Q. By construction Eq is 2-regular, hence consists in a set of disjoint closed 
loops that cover Vd^ (a loop gas; note that not all spanning collections of loops 
on DO are obtained in this way). These loops separate clusters in To from 
clusters in Tq, and it is easily seen that the number of loops is C{E'q) ~ C{Eq) + 

CiEl)-l. AsV-\Eo\ = CiEo)-C{El)+l,wehave2C{Eo) = \VHEo\+CiE'o). 
Hence the FK weight of the configuration Eq may be written as: 

eeEo eeEo 

Consequently, up to normalization, an FK configuration on F (or the dual con- 
figuration on F^^) maps to a loop gas Eq on DO with weight 

e'eE'g 

With w'{e') = (^)' if e' = {{vf),{v'f )); and w'{e') - (^)^ if e' = 
((z;/), («/')) (recall that w{ff)w{vv') = q if (//') = (w')t). 

Parafer mions . 

In the g-Potts case, Nienhuis and Knops ([35]) showed that parafermion cor- 
relators could be computed in terms of the loop gas configuration. Consider Vi 
a vertex of F on the the boundary of the face fi, i — 1,2. One may introduce 
disorder operators fi-k{fi), A^fc(/2), with k e 'Z/qZ ^ fixed. The new config- 
uration space consists of Z/gZ valued 1-forms w on F with duj = kSf^ — kSf-^. 
One may apply the random cluster representation to this state space: if e G E' 
is such that uj{e) = 1, then e Eq with probability Pe, where = w{e). 

The resulting random clusters may not wrap around /i,/2- In this case, one 
also obtains a modified (less canonical) Edwards-Sokal coupling. For instance, 
one may consider a reduced graph the vertices of which are the random clus- 
ters, with edges corresponding to edges of E s.t. uj{e) ^ 1. Then the number 
of assignments compatible with the disorder condition ql'^'^^. One then assigns 
a spin uniformly at random on each cluster, and set a defect line from fi to 
f2 which does not cross any cluster (i.e. is drawn on a dual cluster). There is 
an overall factor q depending on whether one quotients the state space by the 
global Z/gZ symmetry. 
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We now wish to evaluate (cr(v2)^o'(i'i) ''Mfc(/2)/^-fc(/i)), where a{v) e Ug = 

2iT, 2iTi{q-l) 

{1, €•>,...,€ 1 } is the spin variable, and i e Z/gZ, k, £ ^ 0. Plainly, the 
only contributing RC configurations are those in which vi , V2 are connected by a 
cluster and /i, /2 are connected by a dual cluster. Equivalently, in the loop rep- 
resentation, there is a loop L passing through (wi/i) and (u2/2)- Configurations 
containing L will be counted with a phase depending of the isotopy type of L in 
C \ {ill , /i , W2 , /2 } • We now explain how to compute this phase. For definiteness, 
let us fix a defect line 7 from /i to /2. 

Consider the half-loop (resp. L^) of L starting from (wi/i) and ending 
at {^2/2) with vi, V2 on its LHS (resp. RHS). Let us orient both half-loops from 
i^ifi) to {v2f2)- One may also orient all other loops arbitrarily. If all loops have 
an orientation, one may define a height function /i : y<> — J- as follows: if 
v,f are adjacent vertices in V/y, h{f) — h{v) = ±-| according to whether the 

loop crossing vj crosses it from left to right or right to left; then the loops may 
be seen as level sets of h. In the present case, the definition is ambiguous due 
to the confiicting orientations of , L~ . Then h may be seen as an additively 
multivalued function, picking a constant ±7r when circling around (wi/i), (^'2/2); 
this is closely analogous to the notion of magnetic charge we discussed for the 
DGFF. 

Take a path from vi to V2 on the cluster that contains them. It crosses 7 
n times (algebraically: crossing 7 from left to right counts +1, from right to 
left counts -1); and a{v2)(j{vi)~^ = exp( ^'^^" ). Specifically, we can look at the 
leftmost (resp. rightmost) such path, that tracks (resp. L~) on its RHS 
(resp. LHS). Each time one of these crosses 7, it corresponds to crossing 7, 
which contributes to J^dh (the height variation h{f2) — h{fi) evaluated along 
7). Hence nn = dh, and: 

(cr(u2)V(wi)~V/c(/2)Ai-fe(/i))potts = u/l(„i/i)^(^„2/2)exp ( f dh\\ 

\ \ 1 J-r / / o(n) 

where |m| = 1 can be made explicit in terms of 7, arg/iiii, arg/2W2- In the 
case where k or £ generates Z/gZ, it is clear that the LHS depends on k, £ only 
through k£. 

Notice that the expression 



l(i<i/i)«(i<2/2) exp (^2is J d/i)^^ 



7 / / 0(n) 

is defined for any s S R/Z, q > 0, and may thus be used to define parafermionic 
correlators for general FK percolation ([43, 44, 39]). 

5. The 6- vertex model 

The 6- vertex (6V model for short) can be phrased for any finite 4-regular graph 
F = {V,E), such as the square lattice (or some portion of it) or the Kagome 
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lattice. Its free energy per site was evaluated by Lieb, based on a Bethe ansatz 
analysis ([31, 2, 38]). A eonfiguration consists in an assignment of orientation 
("arrow") to each edge, in such a way that every vertex is the origin and the 
endpoint of two oriented edges. Thus there are (2) = 6 possible local configura- 
tions around a given vertex (Figure 2). Graphs with higher (even) valencies may 
be considered: for a 6-regular graph (eg the triangular lattice) , one may consider 
an orientation such that each vertex has 3 ingoing and 3 outgoing edges. This 
defines the (configuration space of) the 20-vertex model (as (g) — 20). Note 
that reverting all orientations preserves the 6V condition. 



Fig 2. 6V configurations 



In the 6V model, one associates a Gibbs weight to each configuration, which is 
a product of local factors, read from the local configuration around each vertex. 
In the planar case (or oriented case), one may order edges cyclically around a 
vertex; then there are four (in-in-out-out) configurations (up to rotation: types 
1-4) and two (in-out-in-out) configurations (types 5-6). Types are also numbered 
so that 1-2,3-4,5-6 are exchanged by reversal of all orientations. 

Let us denote wi(u), . . . ^WQiv) the local Gibbs weights around the vertex v 
and /„ = Iv(yj) G {1>6} the type at v of the configuration w; then the Gibbs 
weight of a; is: 

vev 

Let us denote uj the configuration obtained by reversing all arrows in w. One 
usually requests: 

Wi = W2 = a, ^3 = IV4^ = b, W5 = Wg = c 

so that w{lu) ~ w{ll:). On the square lattice, the rotationally invariant weights 
{wi = • • • = u'4, W5 = We) give the F-model. 

Remark that replacing the weights {wi{v)) with {XyWi{v)) has the effect of 
multiplying weights of all configurations by a constant O-uev -^i"' ^ trivial mod- 
ification. Thus there are essentially two free parameters for the weights at v: an 
anisotropy parameter a/b, and another parameter: 

^ _ + 62 _ 
~ 2ab 

We see the underlying graph F as embedded on an oriented surface E. Then 
there is a dual (embedded) graph F^' , such that each vertex of F''' corresponds 
to a face (assumed contractible) of F, and vice versa; edges of F and F'f are 
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in bijection. Associated to a 6V configuration w on F is a (possibly additively 
multivalued) function h = h{u)) on T\ defined as follows. If /, /' are faces of F 
(vertices of F^) and {vv') is the edge of F between /, /' with orientation given 
by uj, then 

^ h(f) = — if vv' crosses //' from left to right 

h{f') — h{f) = if crosses //' from right to left 

The 6V rule ensures that X]i=o ^ ^ifi) ~ if /o, • ■ • i /s are the four 

faces around v E V, cyclically indexed. This identifies the 6V model as a body- 
centered solid-on-solid (BCSOS) model (see [46]). If S is contractible, h is well 
defined, up to an additive constant. Otherwise, for instance when E is a torus, h 
is additively multivalued, i.e. picks an additive constant (an integer multiple of 
tt) when traced along a non-contractible cycle on S. More precisely, the graph 
1-form J = dh : given by 

./((.//')) = M/')-M/) = -^((/'./)) 

is well-defined and closed: 

dJ{v) = J(e) = 

e^dv 

where f is a vertex of T (a face of F^) and dv are the edges bounding v in Ft, 
oriented counterclockwise. If 7 is a closed loop on F^, J = X)ee7 '^^^'^ depends 
on 7 only through its homology class in iJi(E). 

We may then define several natural order and disorder variables, by analogy 
with free field observables. First we have the current variables J{e^) € ±-| 
(which are sometimes thought of as ±1 spin variables). Then one can define 
natural electric correlators of type: 

(exp(za(/i(/') - him) 

(for definiteness, one needs to fix a path from / to /'), with ±a electric charge 
at/,/'. 

A convenient way to introduce disorder (magnetic) operators is to introduce 
edge defects. An edge e split into two half-edges with outward orientation rep- 
resents a magnetic charge ^ (i.e. the height field picks up ^(27r) when cycling 
counterclockwise around this charge); symmetrically, an edge e split into two 
half-edges with inward orientation represents a magnetic charge — i. 

Let us now restrict to the standard case of the square lattice with periodic 
conditions (^-Z^/L, where L is a sublattice, L ^ Z + tTL where n ^ 1 is a scale 

parameter and r G H = {z : $5(2) > 0} fixed), and isotropic weights, so that we 

2 

may set a = h = 1 , and A = 1 — ^ . We also assume that noncontractible cycles 
have even length. 
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The main conjecture on the fluctuations of the six- vertex model, based on a 
Coulomb gas representation ([34]), is that in the small mesh limit (n — > oo), the 
height field converges to a compactified free field with compactification radius 
^ and coupling constant 

8 . ,c 
Q = — arcsmi - I 

i.e. with action: 




where S = C/rZ is the limiting torus. This is expected to hold in the "regime 
III" phase, i.e. when -l<A<lorO<c<2. 

Relation with FK percolation. 

For clarity we stick to the square lattice with periodic boundary conditions 
(and even periods). There is a close relation between FK percolation and the 
six- vertex model, which is best expressed in terms of the Baxter measure, which 
we now describe. 

Starting from F modelled on the square lattice, we consider its dual graph 
F^ and the diamond graph 0, which is itself locally a (scaled, rotated) square 
lattice. Associated to an FK configuration on F and its dual configuration, we 
have constructed a dense 0{n) configuration (loop gas) on £)<). It is convenient 
to represent its segments as quarter circles: in each face of <), two nonintersecting 
quarter circles connect the four edge midpoints. Starting from the self-dual FK 
weights, the weight of the 0{n) is simply (up to normalization) r/ where £ is 
the number of loops and n = See Figure 3. 




Fig 3. First row: graph T (with periodic boundary conditions); its dual pt (dashed); (}. Second 
row: a primal FK configuration; dual configuration; dense 0{n) configuration 
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The 0{n) state space may be enriched by looking at configurations of oriented 
loops. The winding of each loop if ±27r for contractible loops (which are simple) 
and for non-contractible loops (on the torus). Each quarter circle is either a 
left or right turn and the weight of a oriented loop configuration is set to 

7r 

exp(is — (jright turns] — |left turns])) 

where s is real. Collecting all weights along a given loop L gives e^''^^ "'*"''(^) £ 

{X"\ 1, X}, where x = e^"''- 

Remarkably, this complex measure on oriented loops has two real projections: 
one is a six-vertex model and the other is a weighted FK percolation measure. 

Let us start with the latter: by forgetting orientations, oriented loop config- 
uration configurations map to dense 0{n) configurations, with weight per loop: 
n = x + for contractible loops and n = 1 for non-contractible loops. In the 
thermodynamic limit, one expects the number of these non-contractible loops 
(or, for that matter, the cardinality of any class of macroscopic loops) to stay 
tight. If < (7 < 4, one chooses s so that 2cos(27rs) = ^/g; then the image of 
the Baxter measure projected on dense 0{n) loop configurations has a density 
proportional to loops} | ^^.^ ^-^q random-cluster measure (in 

its loop representation). 

On the other hand, if one retains orientations (at each edge midpoint of <)) 
and forget loops, one obtains a six- vertex configuration on (see Figure 4), 
with weights: 

wi = • • • = CJ4 = 1, tj5 = We = e''"' + e~"^ 
so that a = 6= l, c^ = 2 + y/g. 

\ \ \ \ \ \ 




Fig 4. 6V configurations and corresponding plaguette configurations for oriented loops (on 
each column) 

It is worth pointing out that in this mapping, the parafermionic FK observ- 
ables are mapped to magnetic-electric six-vertex observables. More generally, it 
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is interesting to try to understand which natural observables may be expressed 
in both representations. Let us discuss topological observables ([5, 37, 1]). 

Fix uj = adx + Pdy (thinking of the graph F as having small mesh d and 
being embedded on a macroscopic torus E = C/(Z + rZ), a, j3 fixed, and r G 
i(0, cxd) for simplicity; A,B are the standard cycles on S). Let us modify the 
local Baxter weights as follows: the weight of an oriented quarter circle QT is 
exp(is-|) cxp(i J^^^o;). The weight of an oriented loop is now cxp(±2i7rs) for a 
contractiblc loop and cxp(i u) for a non-contractiblc loop. The possibilities 
for non-contractiblc loops arc: no such loops (in which case either a cluster or 
a dual cluster contains paths homotopic to A and B)] 2k loops homotopic to 
mA + nB, gcd(m, n) = 1 (in which case there are k clusters homotopic to a 
simple path representing mA + nB separated by k dual clusters) . Let us denote 
Pk,m,n the probability of this event under the random cluster measure. In this 
case, summing over orientations yields a relative weight (2cos(77ia + nfi))"^^ . 

On the six-vertex side, each north (resp. west, south, east) pointing edge gets 
a relative weight e*^* (resp. e"'"'', e~*^*, e*""^). Collecting the weights of edges 
along an horizontal (resp. vertical) path on F^^ yields cxp(— J) (resp. 

exp(*^LJ))- 

Summarizing, setting 

/(a,/3) =Po + Pfc,™,„(2cos(ma + n/3)g-i/2)2fc 



2 

A:>1 gcd(m,n) — 1 



we have the exact identity: 



^("'^)-Ee.rexpr2.^^/j_2.^/ J 



/(o,o) V V ^ 



B " JA 



Admitting that the QV height function converges to a compactified free field 
expresses the RHS as a theta function. 



The free fermion point. 



For = -f- 6^, the (square) six- vertex model may be mapped exactly on 
a (weighted) dimer model on its medial graph. Fan and Wu ([10]) first showed 
the existence of a such a dimer representation in the more general context of 
the eight-vertex model. The mapping presented here is in particular apparent 
in the study of the six-vertex model with "domain wall boundary conditions" 
([22]) and its relation with alternating sign matrices and the "aztec diamond" 
problem for dimers ([8, 9, 12]). The correspondence between dimer and 6V 
configurations is represented in Figure 5. Specifically, we start from F which 
is modelled on the square lattice and carries a 6V configuration. Its (planar) 
derived graph DT is also modelled on the square lattice; it has a 2-coloring, 
with black vertices corresponding to horizontal edges of F and white vertices 
corresponding to vertical edges of F. A dimer configuration (or perfect matching) 
on DT is a subset m of edges of DV such that each vertex of DT is an endpoint 
of exactly one edge in m. 
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p ,Q Q a a ,Q 

■-*r\-p-- ■-•t^V- ■-•-f*-- 

b' a TD b' D' a' 

4- 

Fig 5. configurations and corresponding local dimer configuration (bold edge: dimer) 

In the 6V-r-?> dimer correspondence, if {hw) E m, with b (resp. w) the midpoint 
of an horizontal (resp. vertical) edge of F, then: bw points in the direction of the 
edge through b (with its 6V orientation) ±j; points in the opposite direction 
of the edge through w (with its 6V orientation) ztj. It is easy to see that this 
prescription maps dimer tilings to 6V configuration, and that the tiling may be 
reconstructed from the 6V configuration up to local ambiguities (see type 5 in 
Figure 5). One may also introduce defects: for instance, inserting edge defects at 
b, w (with respective magnetic charges —5,5) corresponds exactly to considering 
dimer configurations with monomer defects at b^w ([13]). 

In the dimer model ([27]), each dimer configuration receives a weight 

e£m 

where the We's are fixed edge weights. Here, set We = cos{9) for SW — NE 
edges and We = sin(6') for SE — NW edges. Plainly, the corresponding 6V 
configuration has vertex weights: 

(a : 5 : c) = (cos(6i) : sin(6l) : 1) 

i.e. we recover the free fermion point of the six-vertex model: + b^ = c? or 
A = 0. 

In this case, we see that, up to normalization, the dimer height function (as 
defined in [45, 27] for the dimer model on a bipartite planar graph, taken here on 
every other face) is the height function of the corresponding 6V configurations. 
In this case, building on Kasteleyn's determinantal solution of the dimer model 
([24]), one can verify that the thermodynamic limit is described in some de- 
tails by a (compactified) free field with coupling constant g = 2. In particular, 
one can handle infinite volume Gibbs measures ([3, 28]); current correlations 
and scalar field limit ([26]; for current correlations, it is convenient to consider 
height variations on the even sublattice of F'''); compactified free field limit (for 
an underlying periodic graph, [6]); magnetic and electric correlators ([6]). 
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